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Abstract 

We have investigated the transition from strong to ultra-strong coupling regime between a mid- 
infrared intersubband excitation and the fundamental mode of a metal-dielectric-metal microcavity. 
The ultra-strong coupling regime is demonstrated up to room temperature for a wavelength of 
11.7 fim by using 260 nm thick cavities, which impose an extreme sub- wavelength confinement. By 
varying the doping of our structures we show that the experimental signature of the transition to 
the ultra-strong coupling regime is the opening of a photonic gap in the polariton dispersion. The 
width of this gap depends quadratically on the ratio between the Rabi and intersubband transition 
energies. 
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Intersubband polaritons are quasi-particles issued from the coupling between an intersub- 
band (ISB) excitation of a two-dimensional electron gas and a microcavity photon mode 
In the mid infrared region, ISB polaritons have been demonstrated with planar waveg- 
uides (l-sj, by coupling the ISB transition with a continuum of photonic states. Using these 
systems, electroluminescent devices operating in the strong coupling regime up to room 
temperature have been realized [2[. In the THz frequency range, the strong coupling regime 
has been achieved by coupling the ISB transition with the fundamental mode of a metal 
dielectric metal (MDM) microcavity _ 4]. These cavities are widely used to realize mid and 
far infrared quantum cascade lasers p, lo] . This system allows a highly sub- wavelength con- 
finement of the light as the thickness of the cavity, L, is a small fraction of the wavelength 
(A/L« 1/100) [7]. The coupling energy 2HQr, which is inversely proportional to yL, can 
be therefore strongly enhanced and made comparable with the ISB transition energy E 12 . 
Recently, [8] the record value 2HQr/Ei 2 = 0.48 has been achieved for a transition frequ ency 
of 3 THz. In this case the system enters the so-called ultra-strong coupling regime [9|, 1 1 Of ] . 
characterised by the fact that the anti-resonant and the quadratic terms of the light-matter 
interaction cannot be neglected. An experimental signature of this regime is the opening of 
a photonic gap in the polariton dispersion. In this work, we demonstrate that MDM cavities 
can be used also in the mid infrared wavelength range (5 -15 fim) to achieve the ultra-strong 
coupling regime. By studying three samples, based on identical quantum wells but with a 
different doping level, we have investigated the transition from the strong to the ultra-strong 
coupling regime, up to room temperature. 

Our system is composed of a L = 259 nm thick semiconductor layer, containing 10 
doped GaAs/Al .35Ga .65As quantum wells (9 nm/14 nm), sandwiched between two golden 
plates, one of which is patterned as a strip grating (see inset of figure dH). The MDM 
structure supports photonic modes, whose energy can be tuned by changing the width, s, 
of the metallic strips pjj]. The grating period p, which in our structures is equal to 1.5 /im, 
determines the coupling efficiency with the outside radiation. The energy difference between 
the first two energy levels is £" 12 = 106 meV. In order to study the strength of the coupling, 
we realized three samples (HM2892, HM2951, HM2952) with the same quantum well but 
with different doping levels. Samples HM2892 and HM2951 are delta doped in the middle 
of the barriers, while for HM2952 the Si doping extends in the barriers for a length of 12 nm 
(1 nm of AlGaAs at each side of the quantum wells is left undoped). Shubnikov-de Haas 
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measurements were performed at 1.5 K, resulting in the following values of the electronic 
densities: 8 x 10 11 cm" 2 (HM2892), 1.5 x 10 12 cm- 2 (HM2951), 3 x 10 12 cm" 2 (HM2952). 

The samples were characterized by measuring reflectivity spectra at an angle 9 = 10°. The 
spectra obtained at room temperature for sample HM2952 are shown in Fig. la, for different 
values of the strip width. Two minima are observed in each spectrum, corresponding to the 
upper and lower polariton states. Their energy position varies with s, as this parameter 
determines the detuning between the cavity mode and the ISB transition energy. In the 
figure is also illustrated by a dashed line the transmission spectrum of the quantum well 
measured on an unprocessed sample at room temperature in the multi-pass geometry. A 



minimum is observed, corresponding to the ISB absorption at energy E\2 = y Ef 2 + i? 2 , 
where E p is the plasma energy. The energy £ , 12 corresponds to the limit of the energy of 
the upper polariton for a cavity width infinitely large. It is important to notice that in our 
devices the cavity modes do not depend on the incidence angle 0, as the modes confined 
under the metal strips are decoupled [j^] and give rise to a flat polariton dispersion with 
respect to the component of the photon wavevector perpendicular to the strips, k x . This is 
shown in figure [TJ), where the reflectivity minima measured at 300 K on the same sample 
HM2952 at resonance are shown as a function of k x (white dots). The experimental points 
are compared with the simulated reflectivity, which is presented in a grey scale. In the 
simulation, the quantum well layers are modeled within the random phase approximation 



by a diagonal dielectric tensor, as in Ref. LL3|, while for the gold dielectric constant we used 
the same values as in Ref. UjJ. 

The polariton dispersion extracted from the reflectivity spectra is shown in figure [2^, 
where the minima (bullets) are plotted as function of the cavity mode energy. The minimum 
splitting between the upper and lower branch, 2HQ j r ) is 21.5 meV. Thanks to the very 
high confinement of the light in the MDM cavity, we obtain at room temperature a ratio 
HQr/Ei2 = 0.1, with a system containing only 10 quantum wells. This value is comparable 



to that obtained in Ref. LL4| using a planar microcavity with 70 quantum wells with the same 
doping level as in our sample. The polariton dispersion curve can be calculated by solving 
the following eigenvalue equation [81]: 

(E?-E? 12 )(E?-E%)=4(m R ) 2 El (1) 
where 2HQ.R = ^/J^E p , f w = 0.35 is the overlap factor between the quantum wells and 



the cavity mode and E c is the energy of the cavity mode. Equation [T] provides the eigen- 
values of the complete light-matter interaction Hamiltonian, including the quadratic term, 
which accounts for the polarization self-interaction, and the antiresonant terms, necessary 
to describe a system in the ultra-strong coupling regime [10]. In this case the minimum 
energy splitting 2HQr does not occur at resonance, E c = Ey^ but when the energy of the 



cavity mode is equal to E m [ n = E\ 2 y 1 + ^f- (1 — f w )- Interestingly this value corresponds 
also to that of the asymptote of the lower branch for E c 3> Ei 2 . In the opposite limit, 
for E c <C E12, as already mentioned the upper polariton branch tends to E\ 2 . In Fig. [2] 
we indicated with a star the renormalized ISB absorption energy as obtained from trans- 
mission experiments on the corresponding unprocessed sample (see dashed line in Fig. []]). 
Notice that they are in excellent agreement with the asymptotic value of the polariton for 
E c going towards zero. It is important to remark that the parameters entering in eq. [T] are 
determined from experimental data. The plasma energy is inferred from the value of the 
minimum splitting, while the ISB transition energy is extracted from the value E m [ n of the 
asymptote of the lower polariton dispersion. Moreover, the value of £12 can be consistently 
checked by using the upper polariton horizontal asymptote. The value that we found for E\ 2 
is very close to that calculated using Schrodinger equation in the Hartree approximation. 
From the plasma energy, we extract a value of 8.5 x 10 11 cm -2 for the difference between 
the electronic populations of the two subbands, N\ — N 2 . This value is lower than the one 
expected in our sample at 300 K (2.3 x 10 12 cm -2 ). This discrepancy was already observed 
in the literature [j^] and it can be due to the localization of a fraction of the electrons in 
defects originated from the wafer bonding process. The calculated polariton dispersion, us- 
ing these parameters and eq. dj is shown in Fig. by a continuous line and reproduces very 
well the experimental data including the energy of the measured ISB absorption. Figures [2b 
and [2b presents the reflectivity minima (bullets) and the polaritonic dispersion calculated 
using eq. [1] respectively for sample HM2951 and HM2892. As in fig. [2^, we reported with 
a star the value of the ISB absorption energy extracted from room temperature absorption 
experiments on corresponding unprocessed samples (labeled as ISB abs. in Table Hj). Table [I] 
reports the values of the plasma energy and of the ISB transition energy extracted from po- 
lariton dispersions. In the table we also included the results obtained at 77 K on the three 
samples, where we have observed a small increase of the value of the Rabi splitting, which is 
consistent with an increase of the electronic population on the ground subband. Figure [2^ 
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shows that the polaritonic dispersion presents a photonic gap, delimited by two dashed lines, 
which are the horizontal asymptotes of the polariton branches. As discussed in ref. 8, this 
gap directly results from eq. [Hand has its origins in the quadratic and anti-resonant terms of 
the light-matter interaction Hamiltonian. As a consequence, the observation of a polaritonic 
gap in the dispersion curve can be taken as the experimental signature of the ultra-strong 
coupling regime. The value of the energy gap is: 

E g » f w (E 12 - E 12 ) = (2hn R f I (2E 12 ) (2) 

The right side is a development up to the third order in EpjE\ 2 . From this expression 
it is evident that the ratio E g /E\ 2 has a quadratic dependence on Q R /Ei 2 and therefore 
appears only in the regime of ultra-strong coupling. Fig. [3] presents the values of E g for the 
three samples, both at 300 K (squares) and 77 K (circles), extracted from the polaritonic 
dispersions, and plotted as a function of the measured Rabi splitting. The continuous 
line gives the quadratic dependence on the Rabi energy, by using E\ 2 = 106 meV. The 
left hand side of eq. [2] clearly shows the two quantities responsible for the opening of the 
photonic gap, thus for the entering in the ultra-strong coupling regime. To maximize the 
gap one has therefore to optimize f w and increase the difference between the ISB transition 
energy and its renormalized value due to the collective effect of the two-dimensional electron 
gas. In samples HM2951 and HM2892 the gap is not visible in the polariton dispersion. 
In this case the quadratic and anti-resonant terms can be neglected in the light-matter 
interaction Hamiltonian and the polariton dispersion is well simulated by solving the usual 
strong coupling equation: (E — Ei 2 ) (E — E c ) = {Ml R ) 2 . 

In conclusion, we demonstrated the achievement of the ultrastrong coupling regime be- 
tween a mid-infrared ISB transition (A = 11.7 fim) and the fundamental mode of a MDM 
microcavity. The experimental signature of this regime is the opening of a photonic gap in 
the polariton dispersion. The transition from the strong to the ultra-strong coupling regime 
has been investigated by comparing the dispersions of three identical samples with different 
doping levels. As a perspective, MDM cavities could be employed for realizing active de- 
vices operating in the ultra-strong coupling regime. In this case, thanks to the large value 
of the splitting, a selective electrical injection into the polariton states could be achieved, 
avoiding the excitation of the dark states at the bare ISB transition energy, which dominate 
the transport in existing intersubband polariton devices. 



5 



Acknowledgments 



We thank M. Rosticher for help with the e-beam lithography system. This work has 
been partially supported by the French National Research Agency (ANR) in the frame of 
its Nanotechnology and Nanosystems program P2N, Project No. ANR-09-NANO-007. We 
acknowledge financial support from the ERC grant 'ADEQUATE" . 



[1] D. Dini, R. Kohler, A. Tredicucci, G. Biasiol, L. Sorba, Phys. Rev. Lett. 90, 116401 (2003). 
[2] L. Sapienza, A. Vasanelli, R. Colombelli, C. Ciuti, Y. Chassagneux, C. Manquest, U. Gennser, 

and C. Sirtori, Phys. Rev. Lett. 100, 136806 (2008). 
[3] E. Dupont, H. C. Liu, A. J. SpringThorpe, W. Lai, and M. Extavour, Phys. Rev. B 68, 245320 

(2003). 

[4] Y. Todorov, A. M. Andrews, I. Sagnes, R. Colombelli, P. Klang, G. Strasser, and C. Sirtori. 

Phys. Rev. Lett. 102, 186402 (2009). 
[5] K. Unterrainer, R. Colombelli, C. Gmachl, F. Capasso, H. Y. Hwang, A. M. Sergent, D. L. 

Sivco, and A. Y. Cho, Appl. Phys. Lett. 80, 3060 (2002) 
[6] B. S. Williams, S. Kumar, H. Callebaut, Q. Hu, and J. L. Reno, Appl. Phys. Lett. 83, 2124 

(2003) 

[7] A. P. Hibbins, J. R. Sambles, C. R. Lawrence, J. R. Brown, Phys. Rev. Lett. 92, 143904 
(2004) 

[8] Y. Todorov, A. M. Andrews, R. Colombelli, S. De Liberato, C. Ciuti, P. Klang, G. Strasser, 

and C. Sirtori, Phys. Rev. Lett. 105, 196402 (2010). 
[9] G. Giinter, A. A. Anappara, J. Hees, A. Sell, G. Biasiol, L. Sorba, S. De Liberato, C. Ciuti, 
A. Tredicucci, A. Leitenstorfer and R. Huber, Nature 458, 178 (2009). 
[10] C. Ciuti, G. Bastard, and I. Carusotto, Phys. Rev. B 72, 115303 (2005). 
[11] P. Jouy, Y. Todorov, A. Vasanelli, R. Colombelli, I. Sagnes, and C. Sirtori, Appl. Phys. Lett. 
98, 021105 (2011). 

[12] Y. Todorov, L. Tosetto, J. Teissier, A. M. Andrews, P. Klang, R. Colombelli, I. Sagnes, G. 

Strasser, and C. Sirtori, Opt. Express 18, 13886 (2010). 
[13] L. Wendler and E. Kandler, Phys. Status Solidi B 177, 9 (1993). 



6 



[14] A. A. Anappara, S. De Liberato, A. Tredicucci, C. Ciuti, G. Biasiol, L. Sorba, F. Beltram, 

Phys. Rev. B 79, 201303(R) (2009). 
[15] M. Geiser, C. Walther, G. Scalari, M. Beck, M. Fischer, L. Nevou, and J. Faist, Appl. Phys. 

Lett. 97, 191107 (2010). 



7 



TABLE I: Relevant quantities for the three samples at 300 K (77 K). 





2Ml R (meV) E P (meV) E 12 (meV) Ni 


-N 2 (xlO n cm- 2 ) 


E 12 (meV) 


ISB abs. (meV) 


HM2952 


21.5 (23.5) 36.5 (39.9) 107 (105.2) 


8.5 (10) 


113 (112.5) 


112.5 


HM2951 


14 (16) 23.8 (27.2) 106.3 (108.3) 


3.6 (4.7) 


108.9 (111.7) 


108.9 


HM2892 


9.5 (11.5) 16.1 (19.5) 108.2 (106.8) 


1.7 (2.4) 


109.4 (108.6) 


112 
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FIG. 1: (a) Reflectivity spectra measured for sample HM2952 at 300 K for different values of 
the strip width. The period p is 1.5 /iin. The spectrum in dashed line presents the transmission 
measured on an unprocessed sample in the multi-pass geometry. The inset presents a scheme of 
our devices, together with the relevant parameters of the grating, (b) Reflectivity minima (white 
dots) and simulated reflectivity (grey scale) as functions of k x for s = 1.2 /im; the line represents 
the light cone. 
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FIG. 2: Reflectivity minima at 300 K (bullets), plotted as function of the cavity mode energy. The 
continuous line is the simulated polariton dispersion, following eq. [TJ The dashed horizontal lines 
are the asymptotes of the polariton branches. The star indicates the ISB absorption measured 
from transmission experiments at 300 K on an unprocessed sample. 
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FIG. 3: Energy gap extracted from the measured dispersion of the three samples plotted as a 
function of the Rabi splitting at 300 K (squares) and 77 K (circles) . The continuous line represents 
eq. with E 12 = 106 meV. 
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